Chaotic hysteresis in an adiabatically oscillating double well 
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We consider the motion of a damped particle in a potential 
oscillating slowly between a simple and a double well. The 
system displays hysteresis effects which can be of periodic or 
chaotic type. We explain this behaviour by computing an 
analytic expression of a Poincare map. 

Although hysteresis is a quite familiar and ubiquitous 
phenomenon, no general theory achieves to describe its 
many facets. In condensed matter, hysteresis often ac- 
companies a phase transition, which by nature results 
from the cooperative effect of a large number of degrees 
of freedom. This has recently led some workers to analyze 
it by means of Langevin type ||l[ or master equations ||] 
with infinitely many degrees of freedom, and to propose 
various scaling laws for the area of the hysteresis loop. 
A mean field treatment of this problem |^ reduces it to 
an ordinary differential equation for the order parame- 
ter, with some slowly time dependent external parame- 
ter like the magnetic field. Noise can be incorporated to 
the problem. Similar equations appear naturally to de- 
scribe mechanical or electrical systems, as well as lasers 
Q . Hysteresis effects may appear if the equilibria of the 
dynamical system with static parameter undergo a bifur- 
cation, and scaling laws have been found also in this case 

A paradigmatical example is that of a damped particle 
in a slowly varying potential. If the symmetric potential 
depends on a parameter smoothly interpolating between 
a simple and a double well, the static bifurcation dia- 
gramm looks like the inset of Fig. 2. Imagine now that 
the parameter is oscillating slowly between these extreme 
values. The particle starting close to the initially sta- 
ble origin does not react immediately to the bifurcation. 
Rather, it remains for some time close to the now un- 
stable origin, before falling into one of the newly formed 
minima, which it follows until it merges again with the 
origin. Thus this bifurcation delay, which has been ob- 
served in various physical systems 0| , and was rigorously 
analyzed for the first time by Neishtadt ||] , is responsible 
for metastability leading to hysteresis. 

In this letter we address the following question: for a 
periodic parameter variation, which symmetric potential 
well will the particle fall into during each cycle? In the 
overdamped case, it always chooses the same minimum 
(Fig. la). But at low friction, we found that depending 
on the frequency of the parameter variation, the particle 
may also fall alternatively into the left and right equi- 
librium (Fig. lb), or even go from one minimum to the 



other in a random way (Fig. Ic). We observed the same 
phenomenon in an experimental realization of the sys- 
tem, a pendulum on a rotating table. In this work, we 
show that this surprising behaviour can be understood by 
means of a Poincare map, that we compute explicitly to 
lowest order of the parameter variation. We only outline 
the derivation of this fairly complicated, although essen- 
tially one-dimensional map, postponing rigorous proofs 
to a further publication . 




FIG. 1. Position of the particle (or the pendulum) as a 
function of time, each peak corresponding to the particle 
falling into one of the wells. The sequence of visited wells 
may be periodic, biperiodic or chaotic. 

The equation of motion of a damped particle in a 
slowly varying potential can be written 



q + 27q + $'(g,A(rf)) -0, 



(1) 



where the prime denotes derivation with respect to q and 
< e <IC 1 is the adiabatic parameter. 

Introducing the variables x — {xi,X2) — iq,q), this 
equation can be transformed into the non-autonomous 
first order system 



exi = X2 

ex2 = -2^X2 - $'(a;i, A(r)) 



(2) 
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where the dots indicate now derivation with respect to 
the slow time t = et. 

For fixed A(t) = Aq, the dynamics of (||) are weU- 
known. Around the equihbria xi = g*(Ao), X2 = 0, where 
$'(g*(Ao), Aq) — 0, the hnearization of (||) has eigenval- 
ues a± = -7± - </>", with (/)" = $"((7*(Ao),Ao). 
Hence the fixed point is a saddle if </)" < 0, a stable node 
if < (j)" < 7^ and a stable focus if 0" > 7^. With these 
informations, the phase portrait is easily drawn. 

We will consider potentials of the following type: 
A) is an even analytic function of q such that the ori- 
gin O is hyperbohc when A > 0, a node for A_ (7) < A < 
and a focus for A < A_. For positive A, $ has two min- 
ima ±g*(A) which are nodes when A < A+(7) and focuses 
when A > A+. 

The simplest example is the Ginzhurg-Landau poten- 
tial $((?, A) = -iAg2 + (here A_ = -7^, A+ = 7V2), 
q being the order parameter and A the difference between 
the temperature and its critical value. But there is also 
a simple mechanical system which can be described by 
an equation of this type, namely a plane pendulum on a 
rotating table. In the frame rotating with frequency O, 
it experiences a torque —LMgsinq due to its weight and 
a centrifugal torque JJl^ sing cos g, where L is the dis- 
tance between suspension point P and center of mass G, 
/ the moment of inertia with respect to P, and q the an- 
gle between PG and the vertical. Taking as a time unit 
— {LMg/iy^"^, we obtain the equation of motion 
(1) with $'((7, A) = sing[l - (A-h l)cosg], X = - 1. 

Finally, the function A(t) we consider is periodic with 
period 1, and has exactly one minimum and one maxi- 
mum. In order to analyze the behaviour of the dynam- 
ical system, we want to compute the Poincare map in 
the (a::i, a:;2)-plane during one period to dominant order 
in e. This proceeds essentially by following the motion 
of X along its static equilibria. Let us denote by a'j_{T) 
and a'j_(T) the eigenvalues of the linearization around O 
and q*{X{T)) respectively, and use the notation 

a°'*(ri,T2)= f\eaX-{T)dT, 

Jti 

0°'*(ri,T2)= r lTna%*{T)dT, 
5°'*(ti,T2) = r Re(a^* - a°^*)(T)dT. 

We first study (||) in a neighborhood of the origin. The 
linearized system ex = A(t)x can be solved by construct- 
ing a linear change of variables x — S'(r, e)y transforming 
the equation into ey = B{T,e)y, where B is as simple as 
possible. If A{t) has distinct eigenvalues, one can find 
S regular in a neighborhood of e: = such that B has 
exponentially small off-diagonal terms. Complete diag- 
onalization is possible if we only require S to admit an 
asymptotic expansion in e The full system thus 
becomes 



eyi = a+(T,£)yi -I- &+(yi, ^2; r, e) 

£2/2 = a°(-r, e)2/2 + &-(yi, 2/2; T,e), (3) 

with al{T,e) = alir) + 0{e) and b± = 0{\y\^). We 
assume that a^j_(ro) > and define an adiabatic unstable 
manifold y2 = u{yi; r, e) by the equation 

sii = a°_u + - dyula^yi + b+{yi,u)], 

where the initial condition coincides with the instanta- 
neous unstable manifold at r = tq. The solution, whose 
asymptotic expansion can be computed, may be contin- 
ued to all times such that the origin is a saddle or a node. 
We then carry out the change of variables y2 — + V 
and define in a similar way an adiabatic stable manifold 
yi = v{r];T,e). With yi = v{ri) + (^) becomes 

[a;(T,£)-h/3+(C,?7;T,e)]e 
e?7= [a°(T,£)+/3^(^,77;T,£)]77 (4) 

with I3± = 0{C^ +7]^). 




FIG. 2. Fixed points in the (r, a;i)-plane. Full lines rep- 
resent nodes, dashed lines hyperbolic points, and wavy lines 
focuses. The thin lines are orbits of the system. The inset 
shows the static bifurcation diagramm in the (A,a;i) plane 
(thick lines) with an asymptotic hysteresis cycle (thin lines). 

In the overdamped case, when A(r) always remains 
in the interval (A_,A+), the problem is easy to analyze 
because it can be reduced to a one-dimensional one. As- 
sume that A(r) > for < r < Tb and A(t) < 
for Tb < T < 1. One shows that it is always possible 
to transform (|) into (||), with a°_ + (3- < 0. Thus, 
we conclude from the second equation that ry(r) will go 
to zero exponentially fast, and it is sufficient to con- 
sider the reduced equation on the adiabatic manifold, 
£^ = [a+ -I- /3+(f , 0)]^ = g{CT, e), which undergoes a di- 
rect pitchfork bifurcation at r = ©(e), and an inverse 
pitchfork bifurcation at r = Tb -|- 0(e). 

The most important effect of the adiabaticity of the 
parameter variation is that the bifurcation is delayed. 
The orbit starting with a positive ^ at tq € (rb — 1,0) 
reaches the O(e)-neighborhood of the origin at n = 
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To + ©(ellnel). As long as ^ 
[a° (r) +0{e)]^ so that 



1, 



0{s), we have 



^(r)=e(ri)exp-[«°(ri,T)+0(£)]. 

The smaUness of ^ is thus guaranteed as long as a° < 0. 
For Ti < r < 0, Q!° is decreasing; it has increased again 
to zero only at r = ^'(ti) > 0, which is by definition the 
delayed bifurcation time. Afterwards, the orbit quickly 
jumps on the upper branch q* and follows it adiabatically 
until the inverse bifurcation at Tb, where we can prove 
that C(rb) = Ce^/^. The next delayed bifurcation occurs 
at f = ^(t^). 

It follows that the Poincare map, defined by ^(f + 1) = 
T(^(f)), is a monotonic odd function of ^, such that for 
C > exp-ia°(f,f + 1), one has < cie^/^ < T{^) < 
C2C^^^ ■ This implies that the map has three fixed points, 
the unstable origin and two symmetric stable points. The 
positive solution corresponds to the attractive periodic 
orbit of Fig. la, and, if plotted in the (A,^)-plane, to an 
asymptotic hysteresis cycle (see inset of Fig. 2), whose 
area we are able to prove scales like A{e) ~ ^(0) + ce^l^. 

We now turn to the most complicated case, when the 
amplitude of A(r) is large enough for the origin to be 
a focus for r° < t < r^, and for q* to be a focus for 
T*_ < r < T^. We make a Poincare section at the delay 
time f = ^'(rb — 1), which we assume to be in the interval 
(r*,T;) (Fig. 2). 

The variables C = (■?! ^) ^re defined for r ^ [t° , t° ] and 
for 1^1 < d where d is a constant such that x\ is closer 
to the origin than the focus at q*. We want to compute 
((f + 1) = (Ci,r?i) as a function of C(t) = (Co,?Io), in 
the case where S (0, c?) and ?7o = ^{e). From we 
deduce that ^(r) — d sX t = f(^o) + 0{e), where f (^o) is 
the solution of 



a°(f,f) = -£ln(eo/d). 



(5) 



For most trajectories, f is very close to f , but orbits 
which are exponentially close to the origin at f can be 
appreciably delayed, f (^q) is a decreasing function of 
for Cc Co with f((i) — f and f(Cc) = t+j where 

dexp-ia°(f,r;). 

Next, we construct a particular solution of (|^) x{t), 
wich remains in a neighborhood of the upper branch 
q*{T) and such that ^(jt) = Ce^/"*, r7(Tb) = 0. Writ- 
ing X = x{t) + y, we obtain the equation ey = A{t] e)y + 
b{y,T), where A has eigenvalues a*^{T) + 0{e^/'^). The 
nonlinear term b{y, r) may be decreased to order e by 
the change of variables y — z + x{z), which puts the in- 
stantaneous system into normal form. When t ^ t^, 
the linear part of the equation can be diagonalized. 
Then we have to take care of the turning point rjf., 
where the eigenvalues of A{t) cross and S{t) diverges. 
The problem is solved by carrying out, in the interval 



B into Jordan form, so that one obtains Airy's equation. 
Putting together these steps, we finally obtain 

e(Tb) =e>b)+e"*/^sin(0,/e) 

7?(Tb) = e("*-**)/^sin(0,/e + r), (6) 

where 6* = S*{t^, 1) + 0{e^^'^) is a constant and a* — 
a^Co) + Oie^/^), 0, = MCo) + Oie), with 



a*(6)=«*(^(Co),l) 
</'*(?o) = r(r(Co),T;). 



(7) 



Of course C(''b) depends also on 779 via ?/(f ), but only at 
next -to-leading order. For G Kcjrf], @ is the para- 
metric equation of an exponentially squeezed spiral. 
Proceeding in a similar way around the origin, we find 

C(f+l)=f/C(Tb),with 



u = 



/ sin (^) e-''?/^ sin + 9°^ \ 



(8) 



where 4>° = (/'°(r°,r°) -|- 0(e^/^) is the dynamic phase, 
51 = (5°(0,r°) + 0{e^'^), 6^ = 5°(r^,f) + 0{e^l'^), and 
(5f = (5° + (Jf. The geometric phase shifts Qi can be ex- 
pressed in terms of Airy functions (plus corrections of 
order e^/'^). This transformation acts like a rotation of 
angle <^°/e and an exponential contraction along the sta- 
ble manifold. 

Combining (^ and (^ , we finally obtain the expression 
of the Poincare map = 71(^0, e), r\\ = 72(Co,'7o;£) 
with 



T^ = sin — 



Cei/4 + e"*/-sin(^ 



3("--^*-*?)/-sinf ^ 



sm 



(9) 



r2/3l 



a change of variables which puts 



and T2 = 0{&-^°^l^ 

The expression ^ of Ti is valid for > Cc ■ By sym- 
metry, T\ is an odd function of and for j^ol < Ccj it is 
monotonic as in the overdamped case. All variables ap- 
pearing in (^) are independent of 'yo at lowest order 
in £, except a* and which are given by (Q). 

Since ry is exponentially contracted at each iteration, 
we may replace the two-dimensional invertible Poincare 
map by the non-invertible one-dimensional map ^ 
Ti(^o,0). Its graph is oscillating around Ce^/^ sin((/)°/e) 
with increasing amplitude and frequency as ^ \ ^c- One 
can prove that there is a positive constant [i such that 
when sin(0°/e) > e~'^/^, Ti admits only one positive 
fixed point at ^ ~ e^l^ sin(0°/e), which is stable. In this 
case, we obtain a hysteresis cycle of period 1. Similarly, 
when sin (07£) < -e-'^/^ the map has a stable orbit 
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of period 2, corresponding to the particle falling alterna- 
tively into the left and right well. 

The most interesting situation occurs when | sin(0°/e)| 
< e~'^/^. In this case, we observed numerically a great 
variety of behaviours, including period doubling cascades 
and chaotic orbits. These "chaotic" zones occur at inte- 
ger values of fc = 0°/7re, and the width and height of the 
nth zone are proportional to e~'^^^/'^ This prediction 
has been confirmed numerically. 

The map (^) certainly deserves further study. But the 
most important fact to us is that it explains accurately 
the alternance of periodic, biperiodic and chaotic hystere- 
sis, in accordance with numerical simulations (Fig. 3) as 
well as with the laboratory experiment of the rotating 
pendulum. 
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FIG. 3. Numerically computed bifurcation diagramm of 
the Poincare map for the rotating pendulum. For each value 
of e on the abscissa, we plot the asymptotic behaviour of 
^(f + n), n G IN, for one initial condition. Regions with a 
period-1 and period-2 cycle are separated by small chaotic 
zones, the inset showing an enlargement of the second zone 
from the right. 
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